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^D ' Grating-coupler-induced collective intersubband transitions in a quasi-two- 

t^ ■ 

- ' dimensional electron system are investigated both experimentally and theoret- 

ically. Far-infrared transmission experiments are performed on samples con- 



'^ ' taining a quasi-two-dimensional electron gas quantum-confined in a parabolic 



quantum well. For rectangular shaped grating couplers of different periods we 
observe a strong dependence of the transmission line shape and peak height on 

X" 

?H ' the period of the grating, i.e. on the wave vector transfer from the diffracted 

beams to the collective intersubband resonance. It is shown that the line 
shape transforms with increasing grating period from a Lorentzian into a 
strongly asymmetric line shape. Theoretically, we treat the problem by us- 
ing the transfer-matrix method of local optics and apply the modal-expansion 
method to calculate the influence of the grating. The optically uniaxial quasi- 
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two-dimensional electron gas is described in the long-wavelength limit of the 
random-phase approximation by a local dielectric tensor, which includes size 
quantization effects. Our theory reproduces excellently the experimental line 
shapes. The deformation of the transmission line shapes we explain by the 
occurrence of both types of Wood's anomalies. 
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I. INTRODUCTION 

The optical and electronic properties of electron systems with reduced dimensionality 
like quasi-two-dimensional electron gases (Q2DEG's) formed at, e.g., a modulation doped 
semiconductor heteroj unction have been widely studied in the recent past. Especially the 
collective excitation spectrum of such a system has attracted a lot of attention as it represents 
one of the Q2DEG's most fundamental properties |jl|, but also because of potential device 
applications [Q- [Q. Here, both types of collective charge- density excitations, intrasubband 
as well as intersubband plasmons in the far-infrared (FIR) regime have been intensively 
studied for the last twenty years (see e.g. Ref. P). 

It is known (see e.g. |^) that the intrasubband plasmons always exist for wave vectors 
qy larger than the wave vector of the freely propagating light, gy > u/c, where qy = (g^, qy) 
is the in-plane wave vector (gy = |qy|) of the collective excitation (c: vacuum speed of 
light, u: (angular) frequency), assuming the interfaces of the sample to be parallel to the 
x-y plane. But differently, each intersubband excitation is accompanied by two branches 
of dispersion curves due to the polariton effect 0. One branch is located to the left of 
the light line (gy < uj/c) and thus, is the dispersion relation of a radiative virtual mode, 
and the second appearing to the right of the light line q\\ > u/c describes a non-radiative 
normal mode. Commonly, one calls the Q2D radiative virtual intersubband modes the 
collective intersubband resonances (ISR's), collective intersubband transitions or dimensional 
resonances and the non-radiative normal modes the intersubband plasmons. 

Unfortunately, the investigation of the mode dispersion of the Q2D plasmons co'^^^(gy) is 
not directly accessible in FIR spectroscopy. In general, the Q2D plasmons are accompanied 
by different collective intra- and intersubband transitions due to the intersubband coupling 
(ISC). In this case, the resulting spectrum is of hybrid type. Only for weak ISC nearly 
pure intra- and intersubband plasmons occur. In such a situation, intrasubband plasmons 
can only be excited with fields polarized parallel to the heterointerfaces and which have 
wave vectors q^^ > uj/c, whereas ISR's and intersubband plasmons can only be excited with 



electromagnetic fields having components polarized perpendicular to the interfaces of the 
sample ^,0 and which have wave vectors gy < uj jc and q\\ > u/c, respectively. To solve this 
problem, usually the FIR radiation is coupled to the Q2DEG through a metallic grating of 
period d above the electron system p - [|ri|]. This way, different discrete values of the probe 
wave vector k||„ = {kxn, ky) of the diffracted electromagnetic field parallel to the 2D plane, 
kxn = (^/c) sin 60 + (27r/c/)n; n = 0, ±1, ±2, . . . (60: ray angle of the incident light measured 
from the 2;-axis, which is assumed to be perpendicular to the interfaces of the sample), can 
be excited. Assuming the incident light propagation in the x-z plane and the stripes of the 
grating along the y axis, the incident light couples at the discrete wave vectors q\\ = kxn to 
the collective excitations uj^'^^{q\\) in the Q2DEG, provided the grating is sufficiently close to 
the electron gas. Here, the grating coupler serves as a special optical coupling arrangement: 
(i) to convert the in-plane electromagnetic field component Ex{x,z\uj) of the incident wave 
into the perpendicular component Ez{x, z\u!) of the scattered waves and (ii) induces a wave 
vector transfer k^ — >■ k^n parallel to the interfaces of the system. Moreover, grating couplers 
provide unique advantages over other techniques like the use of prisms |T2|, Brewster angle 



orientation |TB|, or multiple internal reflection (wave guide) geometry |T^: normal incidence 
of the FIR radiation and a relatively large component of the electric field along the direction 
of confinement can be easily achieved. 

The aim of this paper is to investigate both experimentally and theoretically the influ- 
ence of the grating coupler on the coupling efficiency and the resulting transmission line 
shape of the Q2D ISR as a function of the grating coupler period. Both, intrasubband 
plasmons (see e.g. Refs. 0- [|TT| for experiments and Refs. [0- ||19[ for theory) and ISR's 



see e.g. Refs. pO|- ||2^ for earlier experiments and Refs. p^ , |25[| for theory) are intensively 



investigated. Here, we wish to focus on the excitation of the ISR. Unfortunately, the up to 
now developed theories on the grating- coupler-induced excitation of collective Q2D plasmon 
modes is restricted to strong approximations, e.g. optically isotropic media, 2D gratings, 
perfectly conducting gratings, simple half-space geometry etc., which are unsuitable to de- 
scribe correctly the optical response of the Q2DEG synthesized in multilayer systems with 



grating. Moreover, as we will point out below, the influence of the grating coupler on the line 
shape is much more pronounced for the ISR than for the intrasubband plasmon and thus, 
needs to be considered carefully. Because the perpendicular fleld component is zero for the 
zeroth-order beam, the grating-coupler-induced ISR's are nonvertical transitions in k-space, 
i.e. involve wave vector transfers of k^n, ^ = ±1, ±2, . . .. Experiments P^J27| indicate that 



the power absorbed at the ISR frequency is a function of the ratio of the grating period d to 
the intersubband resonance wavelength Xsub- Our theoretical investigations are completely 
universal in the framework of local optics and can be applied to any given optically uniaxial 
multilayered structure with gratings of flnite height. 

II. EXPERIMENT 

The experiments are performed on parabolic quantum well's (PQW's) [^. Here, the 
conduction band edge of an GaAs — Gai^^Al^As quantum well is graded such that it turns 
out to have a parabolic shape in the growth direction. This is achieved by a proper variation 
of the percentage of one of the ingredients (namely the Al-content x). An electron which is 
experiencing this proflle 'sees' a potential which could also result from a homogenous positive 
background such as a doping layer. From Poisson's equation and a back-of-the-envelope 
calculation the potential of this (in the case of a PQW fictitious) positive background with 
density n+ can be related to the growth parameters of the PQW by 

_ £o£.j d Eg _ EqEs 8A 

""+" e2 dz^ ~ e2 W^' ^ ^ 

Here, Es denotes the mean static dielectric constant of Gai-xAl^As forming the PQW, Eq is 
the dielectric permittivity in vacuum, A is the energy height of the parabola from its bottom 
to the edges, the electron charge is -e and W is the width of the grown PQW. Once such 
a structure is remotely doped, the donors release electrons into the well which in turn will 
screen the man-made parabolic potential and form a wide and nearly homogeneous electron 
layer. 



The collective response of such systems has been the subject of many different experi- 
ments and theoretical investigations over the last few years (see e.g. Ref. |2^). The most 
surprising effect that has been observed in the early days of the spectroscopy of PQW's is 
that the electron system absorbs radiation independently from the electron density in the 
PQW only at a single well defined frequency, which can be related to the shape of the bare 
confining potential alone. Self-consistent effects like the renormalization of the subband 
energies in the structures due to electron-electron interactions seems to turn out to be of no 
importance for the optical absorption. In a celebrated paper Brey, Johnson and Halperin 
pO| and Yip [^ generalized Kohn's theorem ||32|, originally derived for the cyclotron reso- 
nance in 3D bulk semiconductor structures. They showed that for the case of a bare perfectly 
parabolic confinement potential long wavelength radiation only couples to the center- of-mass 
(CM) part of the Hamiltonian, leaving the relative coordinates completely untouched. The 
reason for this unique phenomenon is that for a bare harmonic confinement and only for 
this the Hamiltonian can be separated into one part containing only CM coordinates and 
one part containing only relative coordinates. The interaction Hamiltonian for the incident 
light can be shown to couple only to the CM part. In other words, long wavelength radia- 
tion excites a mode in the electron system, which is only connected with the CM motion. 
This collective intersubband mode is sometimes referred to as Kohn's mode or more figu- 
rative ^sloshing mode\ Per construction, it corresponds to the classical plasma oscillations 
of a bulk 3D electron gas of density n^ : Up = (ri+e^/(£:o^s''^e)), where m^ is the effective 
conduction-band-edge mass. If we use Eq. (||) relating n^ to the growth parameters it reads: 

/ 8A \^ 

where Q is the confining frequency of the bare parabolic potential: Vo{z) = ^rrie^'^z^. 

The many-particle picture of this phenomenon |]33|j3^ is that the renormalization of the 
(1-0) subband separation frequency Ai^io = fiio(O) — fiio(^2D_EG); where fiio(^2D_BG) = (^i~ 
£q)/Ti is the subband separation frequency of the PQW with the 2D electron density n2DEG 
(sheet carrier concentration), cancels for gy = with the collective frequency shift A^" = 



ujp^{q\\ = 0) — Qio{n2DEG) of the (1-0) intersubband plasmon a;p°(g||). It is shown [0,0 that 
in such a situation this mode is pinned at uj^'^^q^^ = 0) = fiio(O) = Q independently from 
the electron density. Whereas Kohn's theorem states that the absorption spectrum has only 
one peak at u = Q, the mode spectrum of the freely oscillating Q2DEG consists of all types 
of Q2D plasmons. It should be noted that for gy ^ Kohn's theorem is not valid. 

For our experiments we use different samples that have been cut from one single wafer. 
On top of these samples we fabricate metal gratings of different period d and then compare 
their transmission spectra. Any differences in the spectra obtained for the different samples 
thus can be related to the effect of the different grating couplers. The sample used in 
our experiments is a standard PQW (sample PB48 of Ref. ||29|), schematically drawn in 
Fig.l, in which the Al content x was varied during growth between x = in the center 
of the well and x = 0.3 at its edges. To grade the structure, we used the digital alloy 



technique as described elsewhere p8[. This structure is grown on a semi-insulating GaAs 
substrate and consist of the following layers: layer 1 is the 10 nm thick undoped GaAs 
cap layer (cq), layer 2 is a 200 nm thick undoped GaojAlo^As layer (ci), layer 3 is a 17 
nm thick Si doped (A''^ = 5 x 10^''' cm~^) GaQjAlQ^As layer (^0)5 layer 4 is a 4nm thick Si 
doped {Na = 1.11 X lO^^cm^^) GaQiAlQ_3As layer (di), layer 5 is the 20 nm thick undoped 
GaojAlo^sAs spacer layer (si), layer 6 is the 130 nm thick Gai-xAl^As PQW (L), layer 
7 is the 20 nm thick undoped GaQYAlo^As spacer layer (^2), layer 9 is a 4nm thick Si 
doped {Na = 1.11 X lO^^cm"^) GaQiAlo^As layer (^2), layer 10 is a 200 nm thick undoped 
GaQ^jAlo^As layer (02) and layer 11 is the 500 nm thick undoped GaAs buffer layer (6). The 
substrate is a 500/xm thick GaAs wafer. From the growth parameters we calculate the density 
of the fictitious charge of n+ = 7.4 x 10^^ cm~^ which according to Eq. (0) corresponds to 
an expected resonance energy of about huj^{q\\ = 0) ~ llmeV. A semi-transparent NiCr 
electrode on top of the sample serves as a gate to vary the carrier density in the well and 
alloyed In pellets at its corners provide Ohmic contacts to the electron system. On top of 
the NiCr gate a 50nm thick rectangular Ag grating with stripes along the y axis of width a, 
spacing b between the stripes and period d = a + b is deposited. For the different samples 



investigated here, we used grating periods oi d = 4, 6, 10, 20, 27, 40 and 80 /im, respectively. 
The metaUization for all gratings was chosen to be in such a way that the mark-to-space ratio 
(aspect ratio) is close to t = a/b = 1 (one also introduces the mark fraction / = a/d and the 
open space fraction r = b/d). For comparison, also one sample without any grating coupler 
has been fabricated. To avoid Fabry-Perot type interference effects, the sample substrate is 
wedged by a small angle of approximately 3 degrees. 

Magneto-transport measurements revealed a typical sheet carrier concentration of 
n2DEG = 2 X lO^^cm"^ for all samples. To be able to directly compare the experimen- 
tal results for different samples we made sure that during the actual measurement of the 
FIR transmission for all samples n2DEG was exactly the same by eventually applying a small 
correction bias to the gate electrode [^ . 



Our transmission experiments are performed at low temperatures ( T = 2K) with the 
sample mounted in the center of a superconducting solenoid. Experimentally, we determine 
the relative change in transmission —AT/T = [T(0) — T{n2DEG)]/T{0). T(0) is the trans- 
mission of the sample with the well being completely depleted, T{n2DEG) the one at finite 
carrier densities. The spectra have been taken with a Fourier transform spectrometer and a 
Si composite bolometer was used to detect the transmitted radiation which is guided by a 
ten inches long oversized brass waveguide between the sample and the bolometer. 

In Fig. 2 we depict the results of our experiments, where we plot the relative transmission 
—AT/T for seven samples, investigated as a function of the frequency for zero magnetic field. 
As expected from the dipole selection rules, there is no detectable absorption at the frequency 
of the ISR for the sample without grating. With increasing grating period a peak develops 
at the expected resonance frequency of the ISR. Its line shape and peak height, however, 
strongly depend on the grating used in the experiment. For small grating periods d < Xsub, 
where Xsub = 27rc/(v^^u;) is the wavelength of the FIR radiation at the ISR frequency in 
the GaAs substrate with static dielectric constant Eg, the peak has a Lorentzian shape and 
its height increases with increasing period. For d = Xsub, however, the line shape becomes 
totally destorted indicating some interference effect. For even larger grating periods d > Xsub 
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the relative transmission rapidly decreases. For d = 40/im, the resonance is barely detectable 
and for a grating with period d = 80 fim, no detectable resonance is left. It is also interesting 
to note that the sign of the resonant structure changes when the grating period crosses the 
condition d = Xsub- Whereas for d < Xsub the maximum of —AT/T is positive, for d > Xsub 
the maximum change in transmission becomes negative, i.e. it becomes a minimum. 

III. THEORY 

For the theoretical investigation of the FIR transmission spectra of the PQW structure 
with grating coupler we use our theory, developed in Ref . [^ , which is based on the transfer- 
matrix method of local optics for anisotropic media and apply the modal- expansion method to 
calculate the influence of the grating. Any details of the theory may be found in this paper. 
The combination of both methods results in a generally computationally efficient and stable 
formalism of the optical response of multilayer systems with grating ||36|. As shown in Fig. 
3 the PQW under consideration is modelled by a six-layer system z/ = 1, . . . , 6. Each layer 
is, in general, characterized by its dielectric tensor e„'^^(x, to), where a, (3 = x, y, z, and by its 
thickness di, = \zi, — z^_i\. The layer zi < z < zq contains the rectangular-groove grating of 
height h = di and periodicity d = a -\- b. In the grating region we have for the ffiled stripes 
eQ,i(x, iu) = €^{uj)5af3, where ,^ = a if md < x < md + a and ,^ = 6 if md -\- a < x < {m-\- l)d, 
with m = 0,±1,±2,.... In the experiments samples with metal gratings are used. In 
this case we use the local Drude dielectric function: ea{uj) = 1 — ujp^/uj{uj + i'ja) for the 
metallic stripes, where Up^ is the plasma frequency and 7a is the phenomenological damping 
constant, and we assume for the spacing between the stripes Ef, = 1. The semiconductor 
layers of the quantum- well structure are described by two different types of dielectric tensors: 
(i) £ap{x, u) = e„{iu)6af3 m. Zu < z < Zy-i, if the z/th layer is ffiled by an isotropic medium and 
(ii) ei^^{x,uj) = e^^{uj)5ap is a diagonal tensor if this layer forms the QW and contains the 
Q2DEG which behaves optically uniaxial. This is true in the absence of an external magnetic 
field, where we have e^ = e^ ^ 0, e^ ^ but e^ = e^ = e^ = e^ = ^^ = e^ = 0. 



For the dielectric properties of the background of the GaAs and Gai-xAl^As layers we use 
the so-called es-approximation: 6i^{u!) = Sgu, i- e. we neglect the dynamical properties of 
the optical phonons, but include their screening by the static dielectric constant Esu- This is 
true because the frequencies of the optical phonons are well above the frequency of the ISR 
in the studied PQW's. That layer which contains the Q2DEG we describe by a macroscopic 
local dielectric tensor, which includes the size-quantization effects on the electrons in the 
PQW. The dielectric tensor of the Q2DEG (non-local) is derived in the framework of the 
random-phase approximation (RPA) of the current-response scheme in Ref. [§, wherefrom 
the nonvanishing components in the optical limit {q\\ — > 0) follow in the form: 

UJ[UJ + —) 

4-'N^Mi- %\, ), (4) 



where 

_ 2mJ]io 2 /.v 

ho - — ^ — z^o [^) 

with 

rO-2DEG 

zio = / dz(pl{z)z(poiz). (6) 

Jo 

Herein, ipxiz), K = 0,1,2,..., is the envelope wave function of the PQW and we 

have defined the plasma frequency by uq = [n2DEG^'^ / {.^o^sf^eO'2DEG)Y^'^ i where a2DEG is 

the effective layer thickness of the Q2DEG. Further, t\\ and r^ are the phenomenological 

longitudinal and transverse relaxation times, respectively, VLiq = f^io(^2D_BG) = (^i ~ ^o)/^ 

is the subband separation frequency of the two lowest electric subbands of the effective 

confining potential (non-parabolic) of the PQW and /lo is the oscillator strength of the 

transition 0^1. For the PQW under consideration the subband (bottom) energies Sk , 

the envelope wave functions ifxiz), the Fermi energy Ep and the oscillator strengths Jkk' 

are calculated self-consistently in the framework of the Hartree approximation using the 



method described in our recent paper [^ 
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Assuming monochromatic electric and magnetic fields, E(x, t) = E(x, cu) exp{—iujt) and 
H(x, t) = H(x, a;) exp^—iut), respectively, these fields are given by the wave equations 

V(V ■ E(x,Lu)) - (V ■ V)E(x,cu) = — r(x,w) E(x,cu), (7) 

and 

V(V-H(x,u;))-(V- V)H(x,o^) = -iueoV x [r(x,cu)E(x,c^)]. (8) 

We assume that the incident plane wave travels in the half-space z > zq {u = 0), filled 
by vacuum, within the x — z plane in the negative z-direction with the wave vector k*^") = 
{kx,0, —kf'>), having the components k^ = {u/cjsinQo and k^^^ = {u/cjcosQq. In this case 
the here considered p-polarization (TM-waves) with 

E{x,z\u!) = [E^{x, z\uj),0, Ez{x, z\uj)] (9) 

and 

H(x, z\uj) = [0, Hy{x, z\u), 0] (10) 

is independent from the s-polarization (TE- waves). 

According to the Floquet-Bloch theorem we have in the grating region E*^^)(x+?7i(i, z\uj) = 
exjp(ikxmd)E^^\x, z\ij), where the Bloch wave vector component k^ is defined in the first 
Brillouin zone: — vr/rf < k^ < n/d. Whereas in a multilayer system of homogeneous layers 
k|| = {k^, ky) is a conserved quantity of the whole system, the periodic structure of the 
grating coupler produces an infinite number of propagating waves with ky^ = (k^nyky), 
where kxn = kx + Gn and Gn = ^n, n = 0, ±1, ±2, ... is the reciprocal lattice vector. In the 
presence of the grating the refiected and the transmitted beams are represented by Fourier 
series (Rayleigh expansion [p^^: 

oo 

'Ei{x,z\uj) = ^ exp{ikxnx)Fjn{z,u!), if 2; > ^0 and z < zi. (11) 

This ansatz has to fulfil the wave equation in each layer. In the region z > zq we have 
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oo 

+ Y. Bl^^exp[tiKnX + k^^^z)]i-k^^^e^ + Kne,)}, (12) 

n=— oo 

where A^^^ is the field amphtude of the incident wave, B^ is the amphtude of the nth 
diffraction-order reflected wave and Bq, is the unit vector along the a-axis. From the dis- 
persion relation in vacuum fc^^^ = |k^°^| = uj/c it follows k^^J = [tu'^/c^ — k^^Y^'^. For a layer 
fllled by an isotropic semiconductor we have 



EH( 



Y, exp{ik^nx){A^p,}exp[-ik[':J{z - z^_i)] 






(13) 



where k^'^ = [e^{uj)uj'^/c^ — k^^^^"^ is valid, where A^^^ is the amplitude of the nth order 
diffracted wave propagating downwards and B^^ of that propagating upwards in layer z/, 
and for the anisotropic layer contaning the Q2DEG it follows 



E('^)(x,;z|u;) 



/^O Exx ioj) n=-oo 



Y exp{ikxnx){A^pJexp[-ik'f\z - z^^^i)] 



^\^zn ^x ~r i^xn^z) \ B^^ eXpylK.^^ \Z Zi,-i)\\ K^^ G^ + kxn^z) J 

With fcH = [4^i(^)^7c^ - 4^i^L/4'i)]^/^- 



In the grating region we represent the flelds by the modal- expansion method 38,39 



(14) 



I.e. 



represent the electromagnetic flelds as a sum of the eigensolutions of the wave equation in 
the grating layer. One solves the wave equation in each stripe of the grating and requires the 
boundary conditions [Ez] = and {Hyj = 0, where [A] denotes the change of A evaluated 
at the interface. This results in the dispersion relation of the modes 

{£aPb? + i£bPa? 



-sin{(3aa)sin{(3bb) — cos{f3aa)cos{f3bb) + cos{kxd) = 0. 



(15) 



Equation (|T3p determines for a given pair {uj,kx) a set of eigenvalues {kzi }■ The electro- 
magnetic fleld in the grating layer is the sum over all eigenfunctions: 



E^'^\x,z\uj) 
+X.i{x)(3^i 



c 



(1) 



Y£dX+i{^)kzi 



i(l)„-ifc(^'z 



7Wjk 



■^(i)g-.4^ + EWe^'^l^' 






(1), 



(16) 
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where 



D ,e^Pal^^~^'^) ± P ^-iPal{x-'<nd) 



X+ix) 



md < X < md + a, 



(17) 



JJ^^^iM^-md^a) ^ p^^^-il3„{x-md-a) . md + tt < X < {m + l)d 

The electromagnetic fields, given by Eqs. (|12[)-(|T4|) and (0) have to fulfil the electro- 
magnetic boundary conditions {Exj = and {HyJ = at the interfaces between the different 
layers of the multilayer system. The most profitable method to represent the results is the 
transfer-matrix method. The transfer matrix relates the field amplitudes in one layer to that 
of a different layer. Without the grating coupler the electromagnetic field is characterized 
by two amplitudes in each layer and hence, the transfer matrix is a 2 x 2 matrix. But in the 
presence of the grating coupler the electromagnetic field is represented in each layer by the 
Rayleigh expansion with an infinite number of field amplitudes. It is therefore necessary for 
practical calculations to restrict on a finite number of scattered modes: —Umax < "^ < nmax- 
The field amplitudes of the different diffraction orders of the transmitted and reflected waves 
in layer z/, we arrange in form of column matrices {2nmax + 1 dimensional vectors): 



^M _ ( aM aM aM 



■pi ; 



qM _ t ni'^) dH d(^) 



pn„ 



K-'^pO ) -Op-1, -D. 



pi ; 



■' '■p-rimax 1 ■' '■primax > ' 
■^P-n,nax ' -^pUmax ) ' 



(18) 
(19) 



Note that the incident wave contains only the zeroth order A^il^^^ = {ApQ , 0, . . . , 0). The 
resulting matrix equation, which relates the fleld amplitudes in layer u with that in layer 
z/ + 1 , is given by 



where the 2(2n„ 



P^max 



\ y"'max I 



j-U <-j-V2 

P^max P'^max 



\ ( .1(^+1) \ 



'-j-21 '-j-22 

\ P'^max pf^max j 



^(I'+i) 



P'^max 



\ pii'max I 



1) X 2{2nmax + 1) dimensional transfer matrix is 
Tpn™, Jz^, z^ + 1) 



/ ^11 '-J-V1 \ 

PT^7nax P'^max 



^21 ^22 

\ P'^max pf^max j 



(20) 



(21) 
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The subinatrices are given by 

/ 






I pTlfnax 
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[-'p?lmaiJl(2n.maa;+l) 



(22) 



\ l.^pnmo:cJ(2«maa:+l)l ••• ["'pn^ao: J (2nmaa: + l){2n„„:, +1) ^ 

where i,j = 1,2. These submatrices are (2nmax + 1) x [2nmax + 1) dimensional matrices, 
where we have arranged the elements in the following manner: [^pl^^Jii = ['T'pl^^Jn=on'=o, 



I pTlfjiax 



12 



[X!;? 



prir, 



\n=On'=-l 



' L pn„ 



21 



[X2 



pn„ 



Jn=ln'=0) 



I. primaa:J{2"maa; + l)(2nmaa;+l) 



[^n„aJn=n™.n'=n„„,- The transfer matrix can be written as Tp„_^(z/, z/+l) = [P 



P'^max 



y 



\^pnmaA''^yr^ Dp„^^^(z/+ 1), where Ppn„„;^(i^) is the propagation matrix, which describes the 
propagation of the diffracted waves in layer z/, and Dp„^^^(z/) is the dynamical matrix, which 
depends on the polarization of the waves. The different matrices are derived explicitely 
in Ref. |3^, where the reader can find any details. As shown in Ref. I^G] the propagation 
matrix of a homogeneous layer is given by 

\ 



yii-^nax \ / 



I -p^ 



V 



P'^max 





z/ 







VI 



(23) 



Primax^ ) J 



where the elements of the submatrices read 



\Vln^^S^)]in = exp{ik^J'Jd^)Sin 



(24) 



and 






(25) 



Herein, we have arranged the elements of the submatrices in the following form: 



-pi 

' n=0,Z=n„ 



-pl — Pl P>1 — pl P>1 — P>1 P>1 

' 11 — ' n=0,«=0' ' V2 — ' n=0,Z=-l' ' 21 ~ ' n=-l,«=0'- • • ' ' 12n„maa;+l 

V\„ , i9„ ,1 = 'P3_„ ;_„ . The dynamical matrix of a homogeneous layer is cal- 

^''■max\^^'i'raax ~rJ- "' — ''"max i' — 'i-max "^ »^ ^ 

culated to be 



pll'max V / 



^ Vli (U) P12 



P'^max 



pn„ 



■)22 



ly 



\ 



(26) 



y V'^ z/ P;; ^z/; , 
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where the elements of the submatrices are given by 



Kn^J'')U = PpLaMU = s., 



nl 



(27) 



and 



p;l.mu = -[vii_{u)i, 



nl 



oxx 



Snl- 



(28) 



Equations (|23|) to (^) are vahd for any optically uniaxial medium. The corresponding 
expressions for isotropic media follow if one replaces £^^J{uj) = s^^Ji^) — ^^i^)- 
Defining the transfer matrix of the whole sample by 



- pn, 



{0,N+1] 

lax ^ ' ' 



T?. _(0, 1)T;;„_(1, 2)Tp„_(2, 3) . . . T,„_(iV, iV + 1), 



-pn,, 



(29) 



pTT-max 


\ 


= T^„ (0, N+1) 


y"'max 


\ 


\ '~^pnmax 


J 




\ '~^pnmax 


1 



then field amplitudes {-A-fn^^^i^fn^ax) ^^ ^^^ half-space z > zq are related to that in the 
substrate z < zn hy 



(30) 



Herein, the transfer matrices T£, (0, 1) und T£, (1, 2) result from the two in- 
terfaces of the grating with the homogeneous media. Further, the following is 
valid: The submatrices which result from the interface between two homogeneous 
media are diagonal, whereas those resulting from the interfaces between the grat- 
ing and homogeneous media are non-diagonal. Also the last one can be represented 



in the form T%_{^,\) = [Pp„_(0)]-HD,„_(0)]-iD^c;„_(l) und T^^„_(l,2) = 
[P^„_(l)]-HD^„_(l)]-iD,„_(2). The matrices [Pp„_(0)]-i and [D,„_(0)]-i are the 
inverse of the propagation matrix and of the dynamical matrix of the half-space z > zq 
{u = 0) given by Eqs. p3| ) and (pUf), respectively, and Dp„^^^(2) is the propagation matrix 
of layer z/ = 2 given by Eq. p6|). In the grating region itselfs the propagation matrix 

pG 

prim: 



^^(0) has the same form as given in Eqs. ( p3D to ([25|). But the dynamical matrix in the 



grating region is different from that describing homogeneous layers. The dynamical matrix 
of the grating is given by 
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with 



PH'TYiaX ^ ' 



I^^y' [v^ VIZIS'") 






pl^ma 






(31) 



[^SLMU = [^SLM]n^ = r^^, 



in 



(32) 



and 



[^pn_(^)]' 



nl 



-p?nZs-)u = Ti,k^:;^. 



(33) 



The different element of these matrices are 



1 fd 



a Jo 



{I - exp[ia{f3ai - Kn)]} .^ {1 - exp[-ia{(3ai + Kn)]} , 

^-t^a/ n ; '^^al ~ ~~r ^ 



+ tD^ 



Hal i^xn 

{1 -exp[ib{Pi,i -Kn)]} 



-tF, 



Pal ~r rvxn 

{1 - exp[-ib{pbi + Kn)]} 
exp{ik^na){Pbi + kccn) 



exp{ik^na){f3bi - fc^n) 

a JO e^ 

{l-exp[m(/3ai - A;^„)]} {1 - exp[-ia(/?aj + /c^„)]} 

^-f^a/ 77^ ; ^ "i-Tal 77\ — ; ^ H 



Bailee 



al i^xn) 



Ba[yal + kx 



{1 - exp[ib{Pbi - kxn)]} • i;. {1 - exp[-ib{pbi + k^n)]} 
+ '''^bl 7-, ^ — 77^ ; T i-Tbr 



exp{ikxna)eb{(3bi - k^n) 



exp{ik^ria)eb{Pbi + k^n) 



(34) 



(35) 



Herein, the matrix elements T]^ {n = 0, —1, 1, ...Umax, ^ = 1, ■■■'^nmax + 1) are arranged in 



the form: 1^^ — F^^j^ ,^^q, ^^^2 — ^i=i,n=-ii ^21 — Ti^2,n=oi ■ ■ -^ ^1 



lZ7ltnax~r^ 



/ — 1,71 — TLmax ' ' 



pi = rL 



In the absence of the grating the here presented formalism reduces to the well-known 
2x2 transfer-matrix method, i.e. each submatrix of the transfer matrix given in Eq. (^Of ) 
becomes a single complex function or number and the field amplitudes, Eqs. ([T^) and 
reduce to only one (that of the zeroth order), respectively: 



/A 



pO 






(^+1) 



V-^po 



/ 



\ 



\^pO 



(36) 
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with 



Tpo{v,v+l) = - 



txx rvz txx rvz 



(37) 



C-xx rvz C.XX rvz 

For a multilayer system with a grating coupler, however, transmitted waves of the order 






n = 0, ±1, ±2, . . . ± rimaj: occur. The quantity measured in the experiments is the time- 
averaged power transmission coefficient, which is calculated to be 

^£sN+l\^po\ ^zO n,n'=-n,nax 

where the symbol * means complex conjugate and 

4>nn'{x, Z) = X{kxn' - kxn) + {z - ^iv) (A;i^+'^* - ^-^'^ (39) 

The transmission coefficient depends on both x and z because the grating produces an 
diffraction pattern along the x-direction. It is important to note that not all the orders of 
diffiacted waves can propagate through the sample and from the the sample surface to the 
detector. Only those waves for which k^^^ is real and positive, i.e. if ^ > (^sinOo + ^nY 
are propagating waves above and below the sample. In the case that k^^J is pure imaginary, 
i.e. if ^ < {^sinQo + ^nY the corresponding wave is an evanescent wave with the decay 
length Ln = {2\Imk^^J\) and thus, exists only in the near-field of the grating. Because in the 
experiment the distance between the sample and the detector is usually much larger than 
the decay lengths of the evanescent waves, which cannot transport energy in the negative 
z-direction if the medium below the sample is non-absorbing, only the propagating waves are 
detected. But it is important that the evanescent waves couple with the collective excitations 
of the Q2DEG inside the PQW, which usually is in the near vicinity of the grating, leading 
to a resonant decrease of the reflected and transmitted radiation. Hence, the measured 
transmission coefficient contains via the boundary conditions at the interfaces of the grating 
with the homogeneous layers, which couple the different diffiacted orders, all information 
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about the evanescent higher-order diffracted waves. In the typical situation only the zero- 
order diffracted wave is propagating in the sourrounding. If in this case, the influence of 
the evanescent waves below the sample on the transmission coefficient of Eq. ( p9| ) can be 
neglected. Then, the transmission coefficient reads 

T = \^ ^-^ (40) 



With increasing period d of the grating it becomes possible that the first-order diffracted 
waves become propagating. In the vacuum above and below the sample this takes place 
at (i = Aq, where Aq = 2n/k^^^ = 2ttc/uj is the vacuum wavelength at the frequency of 
the ISR. If this becomes true the transmission spectrum as well as the reflection spectrum 
show the so-called Rayleigh anomaly [p7| , ^0| -|43|l . It consists on a rapid variation in the 
amplitudes of the diffracted orders corresponding to the onset (evanescent -^ propagating) 
or disappearance (propagating -^ evanescent) of a particular diffracted order. This is true 
because with the appearance of a new order of diffracted wave a rearrangement of the field 
amplitudes of the other propagating diffracted waves is caused. The wavelength were this 
takes place is called Rayleigh wavelength X^. But in general, two types of anomalous effects, 
called Wood's anomalies [|^, occur: (i) the above described Rayleigh wavelength type and 
(ii) the resonance-type anomaly. The second type of Wood's anomalies is connected with 



the excitation of a leaky surface wave propagating along the metallic grating ||4l|j43| , |45[| . 

If 2ni + 1 waves are propagating below the sample (zero order plus forward and back- 
ward diffracted waves of higher order), valid ii d > uiXq, it depends on the detector if the 
diffraction pattern is measured or the spatial average. Usually the last one is realized and 
then the transmission coefficient, which only includes the action of the propagating waves 
below the sample, reads 

1 n=ni 

T- - V \AiN+l)\2,{N+l) /.-.N 

^~ I /l(0)|2Z.(0) ^ ' P"" ' ■"" ■ ^ ^ 

^sN+l\^pO I f^zO n=-ni 

The experimental and theoretical investigations performed in this work are focused on 
the question about the coupling efficiency of the grating to excite the ISR of the Q2DEG 



synthesized in the PQW with one occupied subband for a given semiconductor heterostruc- 
ture. It is known (for details see [^) that the couphng efficiency depends on the one hand 
side on the grating parameters (period, height, mark-to-space ratio, material) but on the 
other hand side on the properties of the multilayer system (distance between grating and 
Q2DEG and between the Q2DEG and the lower boundary of the sample, material, etc.). 
Here, we investigate in detail the coupling efficiency of FIR transmission in dependence on 
the grating period d and all the other parameters of the the sample remain fixed. 

IV. NUMERICAL RESULTS AND DISCUSSION 

The numerical calculations are done for the six-layer system depicted in Fig. 3 with the 
following material parameters: layer z/ = 1 is a Ag grating layer with Up^ = 5.69 x lO^^s^^ 
and 7a = 7.596 x lO^^s^^ of height di = h = bOnm; layer i^ = 2 is a GaAs layer of 
thickness d2 = lOnm and £^2 = 12.87 and layer z/ = 3 is a Gai_xAlxAs {x = 0.25) layer 
with Ess = 12.21 and of thickness d^ = 241nm; the fourth layer contains the PQW with the 
bare confining energy hQ = llmeV, in which the Q2DEG is synthesized. For simplicity we 
assume for the background dielectric constant of this layer a homogeneous dielectric constant 
with the parameters of GaAs: £^4 = 12.87 and n2DEG = 2 x 10^^cm~'^, ry = 1 x 10~^"^s, 
r^ = 1 X lO^-'^^s and ^4 = a2DEG = 18.7nm. Layer z/ = 5 is 22Anm thick and consists of 
Gai^xAlxAs and the layer z^ = 6 is a bOOnm thick GaAs layer. In the numerical calculations 
the substrate (region z/ = 7) is filled with GaAs and it is assumed to be of infinite thickness 
in difference to the experimental situation. In the experimental sample the substrate is 
wedged with a mean thickness of approximatly 500/im. In both, theory and experiment this 
is done to avoid the Fabry-Perot resonances in the substrate. 

The numerically calculated relative transmission coefficient —AT/T is plotted in Fig. 4 
for different periods d of the grating and fixed mark-to-space ratio t=l. The wave vector 
transfer from the incident wave to the ISR is for d = Afim: k^^ = 1.57 x lO'^cm"^, k^^ = 
3.14xl0^cm-S...; d = 6^m: k^^ = 1.047 x lO^cm-S k^^ = 2.094 x lO^cm-S .. .; d= 10/xm: 
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ka:, = 6.238 X 10^cm-\ K^ = 1.256 x 10^cm-\ . . .; d = 20/im: K, = 3.14 x 10^cm-\ 
K^ = 6.28 X l(}^cm-\ ...; d = 30/im: k^^ = 2.094 x 10^cm'\ k^^ = 4.188 x 10^cm-\ . . .; 
and for d = 40/im: k^^ = 1.57 x lO^cm'^, k^.^ = 3.14 x lO^cm"^, .... It becomes obvious 
from both experiment (Fig. 2) and theory (Fig. 4) that the efficiency to excite the ISR of 
a Q2DEG inside a given muhilayered-quantum-well system increases with increasing period 
d of the grating coupler. This is true for the given parameters up to rf = lOfim in excellent 
agreement between theory and experiment. It is seen that analogous to the experimental 
result, the theoretically calculated line shapes of the ISR peak show at smaller periods a 
Lorentzian form, which becomes more and more asymmetric with increasing grating period. 
For d = 20(1171 we have the largest magnitude of the maximum in — AT/T, but with an 
asymmetric shape. Increasing further the period transforms the maximum to a minimum 
which vanishes for d ^ oo. Please note that in the case of the asymmetric line shape 
T{n2DEG) > T{0) is valid, i.e. more light is transmitted through the sample with the 
Q2DEG. 

Because it is difficult to define something like the coupling efficiency for a strongly dis- 
turbed line shape, we choose the integrated relative transmission F of the detected lines as 
a measure for it. The result is shown in Fig. 5, where we plot the integrated relative trans- 
mission of the resonant structures for the measured and calculated spectra as a function of 
the grating period. For d < A^u^ we observe a steady increase of F with increasing period d. 
Beyond the vertical dashed line which indicates the condition d = Xgub the integrated rela- 
tive transmission decreases very rapidly. It should be noted, however, that our experimental 
setup only integrates over a finite solid angle given by the distance between the sample and 
the waveguide and by the waveguides aperture. This fact might be important especially for 
propagating higher order waves as we point out below. 

In the case considered here for d < 17.5fim only the zeroth-order diffracted wave is 
a propagating wave in vacuum and in the GaAs substrate. But for a grating period of 
d = 17.5/im at cu = 3 X lO^^s^"*^ (z/ = 159. 3cm^^) the first-order diffracted wave additionally 
becomes propagating in the GaAs substrate and for d = 62.83/im this becomes also true in 
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the region filled by vacuum. In the framework of the e^-approximation the onset of the prop- 
agation of the first-order diffracted wave in the GaAs substrate, which is effectively formed 
by the regions u = 6 and 7, is just aX d = Xsub, where Xgub = 27rc/(^/£^cj) is the wavelength 
of the FIR light in the substrate at the frequency of the ISR. In this case it follows from 
^zi = that E^i = and thus, the electric field of this wave is polarized pure perpendicular 
to the interfaces in the substrate and in all the other layers filled with GaAs. The physical 
situation inside the PQW is very similar but becomes complicated due to the resonance 
structure of the complex dielectric tensor, which may vary strongly in the near vicinity of 
the frequency of the ISR. The onset of the propagation of the first-order diffracted wave in 
the substrate is seen in the calculated relative transmission coefficient. The calculated curves 
show for d = lOfim, d = 20fim and d = 30fim the Rayleigh anomaly at cj = 5.3 x lO^^s"^ 
(z7 = 281.5cm-i), cj = 2.6 X lO^^s"^ {u = 138.1cm-^) and cj = 1.7 x lO^^s"^ {u = 90.3cm'^), 
respectively. It becomes obvious from the theoretical and experimental — AT/T curves that 
just under these conditions the transmission peaks begin to deform. For conditions below 
the onset of the Rayleigh anomaly one measures only the zeroth-order beam. But the peak 
in the relative transmission spectrum results from the evanescent higher-order beams which 
couple with the non-radiative intersubband plasmons. This situation is changed above the 
onset of the Rayleigh anomaly, where the first-order diffracted beam couples with the ra- 
diative ISR. Because the collective intersubband excitations are nearly dispersionless in the 
long- wavelength limit, the resulting ISR peak is a sum of the action of all diffracted beams 
with kxi, kx2, kx3, ■ ■ ■■ Our detailed numerical analysis shows that for the infinite thick GaAs 
substrate only the Rayleigh wavelength type anomaly is responsible for the deformation of 
the peak. The asymmetric line shape is thus caused by the superimposition of the ISR peak 
with the Rayleigh anomaly. For the chosen system the resonance-type anomaly is absent 
in the plotted frequency range, i.e. the comparison of the transmission T(0) and T{n2DEG) 
gives no indication for the excitation of a leaky surface wave carried by the grating. It is no- 
ticeable that the diffracted peaks resulting from the excitation of the intrasubband plasmon 
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a;p*'(g|l) at q\\ = k^i, kx2, kxs, ■ ■ ■ are always symmetrically shaped. 

If one would calculate — AT/T for a sample with GaAs substrate of finite thickness, at 
these frequencies the first-order diffracted wave becomes propagating in the GaAs layers of 
the sample and slightly below these frequencies in the Gai-xAlxAs layers. Possible anomalies 
in this case occurring in the optical spectra due the Rayleigh wavelength effect we call 
internal Rayleigh anomalies. However, in this frequency range the first-order diffracted 
wave is an evanescent wave in the vacuum below and above the sample. To answer the 
question whether these internal Rayleigh anomalies are the cause for the asymmetry in the 
experimentally detected line shape or not, we investigated the multilayer system of Fig. 3, 
but assumed the region z/ = 7 to be filled by vacuum. The relative transmission coefficient 
—AT/T is plotted in Fig. 6 for two different thicknesses d^, assuming here a height of the 
grating oih = lOnm. It is to be seen from Fig. 6 (a) that for dg = 4/im the Lorentzian shaped 
maximum, which appears for d = Anm, deforms with increasing grating period, quite similar 
as shown in Fig. 4, it becomes asymmetrically shaped and for d = 40/im the maximum is 
transformed in a minimum. It is to be seen from the spectra in Fig. 6 (a), calculated for 
d = 20/i?Ti and d = SOiim, that the internal Rayleigh anomalies are absent and thus, they 
have no significant influence on —AT/T. Here, the question arises about the mechanism 
which is responsible for the pronounced asymmetric profiles. In Fig. 6 (b) we plotted the 
corresponding relative transmission spectra for dg = 500nm. It becomes obvious that in this 
case the best coupling efficiency is obtained for d = Gfim. For larger grating periods the peak 
height decreases but with a nearly stable line shape. It should be noted that whereas the 
exponential function of the first-order diffracted wave, appearing in the propagation matrix 
of layer z/ = 6, Eq. (0), is nearly one for small d^ in the plotted frequency range, it varies 
rapidly with the frequency for larger d^. Such a rapid variation may cause Fabry-Perot-like 
resonances of the first-order diffracted wave in the multilayer system. Further, we have 
calculated the total power absorption Ap{n2DEG) = 1 — Tp{n2DEG) — Rp{P'2DEG) for the two 
configurations of Fig. 6 assuming a grating period oi d = 3fim. From the comparison of 
Ap{n2DEG) and Ap{0) (see Fig. 7) it becomes obvious that the Ap{0) spectrum shows for 
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(ie = 4:fim a maximuiii which we attribute to the excitation of a leaky surface wave in the Ag 
grating. This resonance superimposes the ISR resonance to give the asymmetric hne shape 
of the peak. For the smaller layer thickness, d^ = 500nm, the resonance associated with the 
surface wave is absent in the plotted frequency range. Thus, we attribute the asymmetric 
line shape for the finite width substrate due to the resonance-type Wood's anomaly. It 
should be noted that similar line shapes of the ISR resonance were observed in Ref. [^ for 
metal-oxide-semiconductor structures of silicon. In this paper the asymmetric line shape is 
assumed to be caused by the superimposition of the ISR with a non- depolarization- shifted 
ISR. Our rigorous grating-coupler theory, however, shows that the ISR is always accompanied 
by a depolarization shift which cannot be "screened" by the grating. Thus, it seems that 
also in this early experiment Wood's anomalies could be responsible for the asymmetric line 
shape. 

V. SUMMARY 

In summary, we have measured and calculated the optical response of the ISR of a 
quasi two- dimensional electron system in a PQW with different grating couplers on top of 
this structure. Using the transfer-matrix method of local optics and the modal-expansions 
method to calculate the influence of the lamellar grating on the electromagnetic fields, the 
calculated relative transmission describes very well the measured spectrum. It is shown that 
the coupling efficiency of the grating coupler to excite the ISR increases with increasing 
period of the grating up to a certain value where the absorption peak starts to deform from 
a Lorentzian shape to an asymmetric shape. This asymmetric line shape is caused by both 
types of Wood's anomalies, i.e. due to the propagating higher-order diffracted waves in the 
sample beyond the threshold of the Rayleigh anomaly and due the resonance-type anomaly 
(excitation of a leaky surface wave in the grating region). In this case one excites both 
radiative and non-radiative collective intersubband excitations. 
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Fig. 1 Schematic arrangement of the layer structure of the sample containing the PQW 
used in the experiments. 

Fig. 2. Experimentally obtained relative transmission spectra for the ISR in a PQW 
using grating couplers of different grating periods. With increasing grating period 
first the height of the occurring peak becomes larger, then the lines become strongly 
asymmetric and for large grating periods the relative transmission nearly vanishes. 

Fig. 3. Schematic arrangement of the geometry of the multilayer system with grating used 
in the theoretical calculations. 

Fig. 4. Relative transmission coefficient — AT/T of the multilayer system in the near 
vicinity of the ISR of the Q2DEG for different periods d of the grating assuming the 
region z/ = 7 to be filled by GaAs. 

Fig. 5. Integrated relative transmission of the spectra shown in Fig. 2 and 4. As long 
as the grating period d is smaller than the wavelength Xsub of the light corresponding 
to the intersubband transition, the integrated relative transmission increases. For d 
> Xsub the integrated relative transmission rapidly decreases. 

Fig. 6. Relative transmission coefficient —AT/T of the multilayer system in the near 
vicinity of the ISR of the Q2DEG for different periods d of the grating of height 
h = lOnm, assuming the region z/ = 7 to be filled by vacuum: (a) d^ = Afim and (b) 
dg = 500nm. 

Fig. 7. Total absorption coefficient A{n2DEG) and A{0) for the multilayer system of Fig. 
6: (a) (ig = 4/im and (b) rfg = 500nm calculated for a grating period of li = 30/im. 
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